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Abstract. Magnetic fields play an important role in almost all astrophysical phenomena including star formation. But due 
to the difficulty in analytic modeling and observation, magnetic fields are still poorly studied and numerical simulation has 
become a major tool. We have implemented a cosmological magnetohydrodynamics package for Enzo which is an AMR 
hydrodynamics code designed to simulate structure formation. We use the TVD solver developed by S. Li as the base solver. 
In addition, we employ the constrained transport (CT) algorithm as described by D. Balsara. For interpolation magnetic fields 
to fine grids we used a divergence free quadratic reconstruction, also described by Balsara. We present results from several 
test problems including MHD caustics, MHD pancake and galaxy cluster formation with magnetic fields. We also discuss 
possible applications of our AMR MHD code to first star research. 

Keywords: cosmology: theory - magnetohydrodynamics - methods: numerical 
PACS: 95.30.Qd 
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INTRODUCTION 

Adaptive mesh refinement(AMR) cosmological hydro- 
dynamics simulations play an important role in the study 
of structure formation of different scales from galaxy 
clusters to first stars in the past ten years. Its ability 
to achieve very high resolution in large scale simula- 
tions with relatively small computer resources has helped 
us to understand the first stars formed in our Universe. 
The possible effects of magnetic fields have been largely 
ignored. It is well established that magnetic fields are 
present on different scales, from intracluster medium to 
interstellar medium. The origin and evolution of these 
magnetic fields and their role on the structure forma- 
tion are still unclear. So, it is imperative to include mag- 
netic fields into current hydrodynamics AMR cosmol- 
ogy code. In this paper, we present the newly developed 
MHD version of the Enzo, which is wildly used in the 
study of first stars IfTl IT4l IT31 . 



MHD WITH ENZO 

The MHD equations in the comoving coordinates are: 
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where all variables have their usual meaning, a is the ex- 
pansion parameter. To track the pressure more accurately 
in the supersonic region, we have also implemented the 
modified entropy equation given in Ryu et al. |Qj5| an d 
the internal energy equation given in Bryan et al. |@] in 
our code. 

The MHD solver used for all the test problems here is a 
high-order Godunov-type finite-volume numerical solver 
developed by S.T. Li Q E). This solver was recently 
successfully used to study magnetic jet problems l9lfT2l 

na. 

We used a constrained transport(CT) scheme flux CT 
O to maintain divergence-free magnetic fields. For the 
AMR hierarchy in Enzo, we used a modified divergence- 
free reconstruction scheme original proposed by Balsara 
including second order accurate divergence-free re- 
striction and prolongation for magnetic fields. Details of 
the CT and AMR of magnetic fields in the MHD Enzo 
can be found in Collins et al. [5|. 

The MHD module has been tested and shown to 
be compatible with other physics packages installed in 
Enzo, such as radiative cooling, star formation and feed- 
back. 



TESTS 



Cluster Formation with Magnetic Fields 



MHD Caustic and Pancake 

The MHD Caustic test is taken from Li et al. [10| 
which generalizes the HD test of lfl6ll . The initial sinu- 
soidal velocity field in the x-direction has the peak value 
1/27C, the initial density and pressure has been set to be 
uniform with p = 1 and p = 10~ 10 . Then caustics will be 
formed because of the compression by the velocity field. 
An initial uniform magnetic field of 10~ 3 in code units 
in the y direction was added to the simulation. Figure [T] 
compares the density and B y at t = 3 of unigrid and AMR 
runs. AMR run had 256 cells in the root grid and 2 level 
refinements by 2. The AMR solution is indistinguishable 
from a uniform grid solution with 1024 cells. 

Pancake is another standard test problem of cosmo- 
logical hydrodynamics simulation [16]. We have run the 
collapse of a one-dimension pancake in a purely bary- 
onic universe with Q. = 1 and h=\. Initially, at a, ■ = 1, 
which corresponds to z, = 20 in this test, a sinusoidal 
velocity field with the peak value 0.65/(1 +zi) in the 
normalized unit has been imposed in a box with the co- 
moving size 64-hr Mpc, so that the shock forms at z = 1. 
The initial baryonic density and pressure have been set 
to be uniform with p = 1 and p = 6.2 x 10~ 8 in the nor- 
malized code units. We applied initial uniform magnetic 
fields By = 2.0 x 10~ 5 G, B x = B Z = in the simulation. 
We did the calculations with unigrid with 1024 cells and 
AMR run with 256 cells of root grid and 2 level refine- 
ments. Figure [2] shows the density and B y at z = 0. 





Cluster formation is one of the problems most widely 
studied by Enzo [ 1 1 1 . We have done this problem and 
compared with results from Enzo-ppm to test our new 
code in large scale structure formation. The simulation 
is a ACDM model with parameters h = 0.7, Cl m = 0.3, 
£2/, = 0.026, £2a = 0.7, a& = 0.928. The survey volume is 
256 hT x Mpc on a side. The simulations were computed 
from a 128 3 root grid and 2 level nested static grids in the 
Lagrangian region where the cluster forms which gives 
an effective root grid resolution of 512 3 cells (0.5 h~ l 
Mpc) and dark matter particles of mass 1.49 x 1O 1O M . 
AMR is allowed only in the region where the galaxy clus- 
ter forms, with a total of 8 levels of refinement beyond 
the root grid, for a maximum spatial resolution of 7.8125 
h~ 1 kpc. 

We first present the results of no initial magnetic 
fields and compare them with the results from Enzo- 
ppm. The cluster parameters from the MHD code is al- 
most identical to those from Enzo-ppm. The virial radii 
are 2.229Mpc from hydro and 2.226Mpc from MHD 
while the virial masses are 1 .265 x 10 15 M Q for hydro and 
1 .260 x 10 15 M o for MHD. Figure [^compare the projec- 
tions of the baryon density and temperature. 
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FIGURE 1. Plots of density and y component of magnetic 
fields of MHD caustics at t=3. The initial magnetic field is 
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FIGURE 2. Plots of density and y component of magnetic 
fields of MHD pancake. The initial magnetic field is 2 x 10 G. 



FIGURE 3. Logarithmic projected gas density(top) and log- 
arithmic projected X-ray weighted temperature(down) at z = 
of adiabatic simulations. The images cover the inner 4 A~'Mpc 
of cluster centers. The left panels show results from the PPM 
solver and the right panels show results from the MHD solver. 
The units of density and temperature are M^/Mpc and Kelvin 
respectively. 

We have performed simulation with initial magnetic 
fields, B x = B z = 0, By = 1 .0~ 9 G with radiative cool- 
ing, star formation and star formation feedback. Figure 
Hlshows the images of baryon density, temperature, mag- 
netic energy density and Faraday rotation measurement 
of the cluster center. The rotation measurement is inte- 
gral along the projection direction. 
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FIGURE 4. Images of the baryon density (logarithmic, 
Mq/Mpc^), temperature (logarithmic, Kelvin), magnetic en- 
ergy density (logarithmic, erg /cm ) and rotation measurement 
(rad/m 2 ) of the inner l/i~'Mpc of cluster center. The initial 
magnetic field is By = 1 .0 G at z = 0. 



Another simulation we performed is without initial 
magnetic field but with the Biermann battery effect. The 
induction equation is modified by adding an additional 
battery term| 6 1 : 



DISCUSSION 

We have introduced our MHD module for Enzo and pre- 
sented some test results. With this new module, we have 
the ability to perform AMR MHD cosmology simula- 
tions. This module uses high accuracy TVD MHD solver 
and employs CT and AMR divergence-free magnetic 
fields reconstruction scheme to guarantee divergence- 
free of magnetic fields. MHD simulations using this 
MHD module can use all the exist physics packages in 
Enzo. 

It is widely believed that the magnetic fields should 
have little effects in the formation of the first stars, since 
maybe there were no magnetic fields at all. But even 
there are no magnetic fields from the early universe, 
the Biermann battery effect should have generated some 
very small fields during the collapse to form those stars. 
If the first generation stars spin very fast, these small 
seed fields could be maintained and amplified by dynamo 
effect. Then the magnetic fields in the first stars may 
pollute the environment by their explosion and play a 
part in the formation of next generation stars. This would 
be investigated in the near future. 
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dt n e e 



(7) 
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where c is speed of light, p e is the pressure of electron, n e 
is the electron number density, e is the electron charge, 
rriH is the hydrogen mass and % is the ionization fraction. 
We took x — 1 constant in space in our simulation. 
We performed this computation with radiative cooling. 
Figure [5] shows the projection of logarithmic baryon 
density and the magnetic energy density of the cluster 
center. 
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FIGURE 5. Images of the baryon density (logarith- 
mic, M^j/Mpc^) and magnetic energy density (logarithmic, 
erg J cm ) of the simulation with Biermann battery effect of the 
inner lfr^'Mpc of cluster center. 
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